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We present a complete analysis on energy and wavefunction of Harmonic oscillator 
with simultaneous non-hermitian transformation of co-ordinate [x —>- P+' A D ) an d 


momentum (p —» 


(jp+ifix) 


y/ (1+/3A) ' 

^/(i +/ 3 A) ; us * n & perturbation theory under iso-spectral condition. 
Further we notice that two different frequency of oscillation (iy 1; tc 2 )correspond to 
same energy eigenvalue ,which can also be verified using Lie algebraic approach [ 
Zhang et.al J.Math.Phys 56 ,072103 (2015) ]. Interestingly wave function analysis 
using similarity transformation[F.M. Fernandez, Int. J. Theo. Phys. (2015) (in Press)] 
refers to a very special case. 
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I. INTRODUCTION 

In quantum mechanics, Harmonic oscillator plays a major role in understanding 
limitations of various approximation methods such as variational method, W.K.B. 
method, perturbation method etc. The basic advantage of this oscillator is that 
its energy levels and eigenfunctions are exactly known [1], However, when the co¬ 
ordinate and momentum simultaneously under goes a non-Hermitian transformation 






[2-4], its energy eigenvalue can be iso-spectral to original Harmonic oscillator. Un¬ 
der iso-spectral condition, the wavefnnctions of the transformed Hamiltonian differ 
drastically from that of the original Hamiltonian. However, a complete picture on 
wavefunction is still a missing subject. Further, no such explicit calculations on wave- 
functions are available at present. Hence the aim of this paper is to present a complete 
picture on wavefunction and energy under iso-spectral condition using perturbation 
theory. 

II. Energy levels and Wavefunction of Simple Harmonic Oscillator (SHO) 

Let’s consider the Hamiltonian of SHO [1] 


■n 2 'r 2 


(i) 


whose exact energy eigenvalues is 


E n — (w + —) 


and corresponding wavefunction is 

r 

= 


-H n (x)e "2 


( 2 ) 


( 3 ) 


y / 7r2 n n 

where H n (x) is the Hermite polynomeal. 

III. SHO under non-Hermitian transfromation of co-ordinate (x) and 
momentum ( p ) 


Let us consider the non-Hermitian transformation of x and p as [2-4] 

x + iXp 


and 


x —> 


p y 


y^T+T^A) 

p + if3x 


( 4 ) 


( 5 ) 


\/(l + Pty 

In this transformation we noticed that the transformed co-ordinate and momen¬ 


tum preserve the commutation relation i.e. 

[x,p\ = i 


( 6 ) 
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Now the new Hamiltonian with transformed x and p becomes non-Hermitian in 


nature and is 


(p + i(3x) 2 (x + i\p) 2 

H = + 


( 7 ) 


2(1 + A/1) 2(1 + X/3) 

IV. Second Quantization and Hamiltonian 

In order to solve the above Hamiltonian (Eqn. (7)), we use the second quantization 
formalism as 

(a + a + ) 


x = 


y/2uo 


and 


P = HI 7,/ a 


( 8 ) 


( 9 ) 


where the creation operator, a + and annhilation operator a satisfy the commuta¬ 
tion relation 

[a, a + ] = 1 (10) 

and u is an unknown parameter. The Hamiltonian in Eqn. (7) can be written as 


h = h d + h n 


( 11 ) 


where 


H d = [(1 - X 2 )uj + 


(1 - (3 2 ) J2ci + a + 1) 
u> J 4(1 + X/3) 


and 

a 2 ( Q+ ) 2 

Hn = U 4(1 +X/3) + V 4(1 +X/3) 

V = [-w(l - A 2 ) + (1 ~^ } - 2(A + (3)\ 

OJ 

U = [-w(l - A 2 ) + (1 ~ ^ + 2(A + (3)} 

UJ 

V.(a). Zero Energy Correction Method (Case Study for U=0) 

Now we solve the the eigenvalue relation: 


( 12 ) 


(13) 


(14) 


(15) 
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H^ n (x) = e n x V n (x) (16) 

using perturbation theory as follows. Here we express 

<* = 4 0) + E 4 n) (i7) 

m= 1 

The zeroth order energy satisfies the the following eigenvalue relation 

H D \tp n >= H D \n>= e^\n> (18) 

where is the zeroth order wave function and e ^ is the mth order perturbation 


correction. 


(o) = ( 2n + 1 ) rn _ \ 2 )u + (1—E>] 
n 4(1 + A/3) 1 J u 1 


E 4"° = + 4 + 4 } +. (20) 

m= 1 

The energy correction terms will give zero contribution if the parameter is deter¬ 
mined from non-diagonal terms of Hn [2] 

Let the coefficient of a 2 is zero [5] i.e. 

U = [-w(l - A 2 ) + (1 ~^ + 2(A + /3)] = 0 (21) 

id 


which leads to( considering positive sign) 


CU = U>1 = 


(1 + P) 

(1-A) 


In this case, 


4"’ = (n + 


Now the perturbation correction term is 


TT = y ( Q+ ) 2 = ( A + P) ( +) 2 

N V 4(l + A/3) 1 + A/3 ’ 
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In this case one will notice that 


< nlH^ln — 2 >= 


Jn(n-1) 
4(1 + A/3) 


(25) 


< n — 2\H N \n >= 0 (26) 

Hence it is easy to note that all orders of energy corrections will be zero. Let us 
consider explicitly corrections up to third order using standard perturbation series 
given in literature[1,2,6-8], which can be written as 

4 1} =< ^n\H N \^ n >= 0 (27) 


e ( 2 ) = y 

n / -j 
k^n 


< 'tp n \H N \'lpk >< Ipis | H N | r !p n > < 'l/j n | H N | lpn+2 >< ^n+21 #Jv| > 




- &) 


(3) _ ^ ^Pn\HN\lpp X lJjp\H]y\'ljjq ><C 'Ipq \ H]\f \ l[)n 5” _ 

n / . 


= 0 


p,q 


or 


(ef-e^e^-e^) 

( 3 ) _ < ^n\H N \^ n+2 >< ^ n +2\H N \^ n+A >< ^ n +A\H N \^ n > _ 


= 0 


= 0 


(28) 


(29) 


(30) 


(JO) JO) \(J 0) (0) \ 

( e n ~ t n+2 ){tn — t n +4j 

Here second order correction is zero due to < '0 n |i3jv|'0n+2 >= $n,n +4 and third 
order correction is zero due to < ^w-kI-^vIVTi >= <S n +4,n+2- Similarly one can notice 
all correction terms el™'- 1 will be zero. Hence 


£n = 4" 1 = =(n+~) 


(31) 


l^n >= ( 


_ / Vgr 


^/^T2 n n\ 




(32) 


with 


< 'Ipnl'tpn > = 1 


(33) 
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Form Eq.(2) and Eq.(23), one can see that the total energy of the Harmonic 
oscillator and the Harmonic oscillator with non-hermitian transformation remains 
the same. 

V.(b). Coresponding Wavefunction using Perturbation Theory 

Here we find the wavefunction as 


J (n + 2)! J (n + 4)! J (n + 6)! 

= I i’n > +/a,/ 3-- 7=7-|VVf2 > +(f\,p) -_ 7=F-IVW-4 > + (/a,/?) -„„ rr l^n+6 > +• 




8 VnJ. 


48 VnJ. 


where ~ com P ac t form one can write, 


*!?> = Z 


+ &) ifc. /( n + 2 ^)' 


k~o ( X + A ^)' 


m 


| , 0n+2/c 


UJl 


(34) 


(35) 


The normalization condition here can be written as[6] 


< </>» i«f > >= i 


(36) 


So also the eigenvalue relation 

<V-„I ff|f<f>>=£„=(n+i) (37) 

VI.(a). Zero Energy Correction Method (Case Study for V=0) 

Let the coefficient of (a + )2 is zero [2] i.e. 

V = [—oj(1 - A 2 ) + (1 -^ - 2(A + P)\ = 0 (38) 

(d 


which leads to 


UJ = U) 2 = 


( 1 -/ 3 ) 


(14- A) (39 > 

In this case, we would like to state that u calculated using similarity transformation 

[3] remains the same as U 2 - Now the perturbation term becomes 


H n = U 


In this case one will notice that 


a 2 _ (A + (5) 
4(1 + A/3) “ 1 +A/3' 


< <t>n\HN\4>n+2 >~ U 


\J[( n + l)( n + 2)] 

4(1 + A 0) 


(40) 


(41) 
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< 4 > n+2\HN\4>n >~ 0 


( 42 ) 


Hence we have 


4 1 ' =< <t> n \H N \4>„ >= 0 


(43) 


4 2) = E 

k^n 


< 4>n\HN\4>k >< 4>k\HN\4>n > < 0n|-^iv|0n-2 >< 4>n-2\H jv|0 n > 


(M _ (oh 
\^ n ^k ) 


(4" 1 - 4b) 


(3) _ 4 4^n | H]\f | ( f>p b* 4 (f)p | H]\f | (f)q > (j)q | HjV 1 0n ^ _ 

n / . ~ 


P,<? 


(ei 0) -6i 0) )(ei 0) -6( 0) ) 


= 0 


= 0 

(44) 

(45) 


or 


(46) 


( 3 ) ^ ^n\H 2 0n—2110n—4 bb *fin—4 1 -^TV | q 

e?l = Tlo) (0) w (0) (0) : = u 

(tn t n _2)\tn t n -4) 

Here second order correction is zero due to < 0 n |Hjv|</> n _2 >= S n , n -2 and third 
order correction is zero due to < 0 n _4|Hjv|0 n >= d n _ 2 , n - 4 . Similarly one can notice 
all correction terms e)™- 1 will be zero. Hence 


e« = 4°’ = En> = («+-) 


(47) 


which is the same as the energy level of harmonic oscillator as given in Eq(2) and 

(48) 


\4>n >= )itf„(v£«:)e-*£ 


' y/ir2 n n\ 

VI. (b). Coresponding Wavefunction using Perturbation Theory 


Here we consider the wavefunction as 


~ \4>n > +f\,p - 7= = I2 > +(/a,/s) 2 - 7===\ ( j ) n-4 > + (/a„s) 3 


VnJ 


V< n - 2 > ! 


8\/( n “ 4 ) ! 


48^/(n-6)! 

(49) 


|0n-6 > +•• 
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In its compact form, one can write 



*>(k) _X^( (^ + ^) \k ij, 

(! + A/3) \J(n — 2k)\2 k k\ 

(50) 

Here we 

notice that for x —> 00 i.e. 



4> n (x —> 00 ) —> 0 

(51) 

and 

$^(o; —> 00 ) —> 0 

(52) 

In this case, the normalization condition can be written as[6] 



< >= 1 

(53) 


So also eigenvalue relation 

<<t>n\H\^)>=E n = {n+ 1 -) (54) 


VII. Comparision with Similarity Transformation using Lie-algebra [4]. 

In the above we notice that two different frequency( wi, w 2 )corresponds to same 


energy eigenvalue .Now we compare our results with that of Zhang et.al [4] using 
Lie-algebra as follows.Previous authors consider the Hamiltonian 

H = So(q~^o T ~) + Si(u~^) J T s 2 (i T .sqnP T s^ci (55) 

having energy eigenvalue 


e n 



- Asis 2 (n 



+ 


S2S3 + S1S4 — S0S3S4 

s 0 — 4 sqSiS2 


(56) 


Case-I w — w\ — 
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In this case we have the Hamiltonian 


tt (l + 1 | 1 \ & + P) / +\2 

H -(aa+ 2 ) 2(1 + A/?) ° 

(57) 

Now comparing we get 


«o = 1 

(58) 

_ (-^ + /5) 

' Si “ 2(1 + A/3) 

(59) 

S 2 = 0 

(60) 

s 3 = 0 

(61) 

S 4 = 0 

(62) 

Hence the 


en = [n + -) 

(63) 

which is same as given in Eq(31). 


Case-II w = w 2 = • 


In this case we have the Hamiltonian 


tj t + 1 \ (A + /3) 2 

H — (tz (xH - ) + . cl 

y 2 2(1 +A/3) 

(64) 

Now comparing we get 


So = 1 

(65) 

(A + /3) 

52 " 2(1 +A/3) 

( 66 ) 

Si = 0 

(67) 

s 3 = 0 

( 68 ) 

S 4 = 0 

(69) 

Hence the 


/ 1 \ 

Cn = (n + -) 

(70) 
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which is same as given in Eq(46). Now we see results of Lie-algebra matches with 
that of perturbation theory on energy level calculation. 


VIII. Comparision with Similarity Transformation [3] 


It is worth to mention that Fernandez [3] has calculated groundstate wave func¬ 
tion of this oscillator using similarity transformation.The explict expression for wave 
function [3] is 


|<£o >= 


v—) 4 e 

7r 


(71) 


It is easy check that our result in Eq. (45) with n — 0 remains the same as that of 
Fermamdez [3]. However, we do not have literature for further comparison. 


IX. Conclusion In this paper, we suggest a simpler procedure for calculating 


energy levels and wave function of the non-Hermitian harmonic oscillator under si¬ 
multaneous transfromation of co-ordinate and momentum using perturbation theory 
. Further energy eigenvalue calculated using perturbation theory matches nicely with 
that of Lie algebric method [4.At this point it is necessary to mention that the ground 
state wave function calculated by Fernandez [3] refers to zeroth order wave function 
as reflected in VI(b).However we present a complete picture on wave function. In all 
the cases, we show that the energy levels remain the same as that of simple Harmonic 

oscillator. Here we would like to state that if the parameter c o is determined using 

, ^ 0 ^ 

variational principle [9] i.e. = 0, then one has to calculate all orders of perturba¬ 
tion corrections i.e. J2m =2 e n^ because < n\H]^\n+2 >^ 0 so also < n+2|iAv|n > ^ 0 
and it will be a cubersome process. However, if the parameter is determined either 
using < n\H N \n + 2 > or < n + 2\H N \n > as done in above procedure then the energy 
levels and wave ftnction are obtained easily. Further we notice that 

< (f) n | Hu | <j) n >=< r tp n | H d | ijj n >=< (j>n\H\^ >=< ^ n \H\^ >= (n+ 1 ) ( 72 ) 

This directly reflects that non-commuting operators i.e. [H, Ho] ^ 0 corresponds 
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to different wave functions [1]. Lastly, one can use nonlinear perturbation series 
[10] and easily get convinced that the parameter c o determined using the condition 
< n\H N \n + 2 > or < n + 2\H N \n >, the calculations become very easy to get the 
desired result. 
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